Heat Transfer between Graphene and Amorphous SiOs 
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We study the heat transfer between graphene and amorphous Si02. We include both the heat 
transfer from the area of real contact, and between the surfaces in the non-contact region. We 
consider the radiative heat transfer associated with the evanescent electromagnetic waves which 
exist outside of all bodies, and the heat transfer by the gas in the non-contact region. We find 
that the dominant contribution to the heat transfer result from the area of real contact, and the 
calculated value of the heat transfer coefficient is in good agreement with the value deduced from 
experimental data. 



Graphene, the recently isolated 2-dimensional (2D) 
carbon material with unique properties due to its linear 
electronic dispersion, is being actively explored for elec- 
tronic applications Ij. Important properties are the high 
mobilities reported especially in suspended graphene, the 
fact that graphene is the ultimately thin material, the 
stability of the carbon-carbon bond in graphene, the 
ability to induce a bandgap by electron confinement in 
graphene nanoribbons, and its planar nature, which al- 
lows established pattering and etching techniques to be 
applied. 

Recently it has been found that the heat generation in 
graphene field-effect transistors can result in high tem- 
perature and device failure[2|- Thus, it is important 
to understand the the mechanisms which influence the 
heat flow. Because of surface roughness the graphene 
will only make partial contact with the Si02 substrate, 
which will reduce the heat transfer coefficient as com- 
pared to the perfect contact case. In Ref. [2| the temper- 
ature profile in the graphene under current was obtained 
by studying a Raman active phonon band of graphene 
(the position of the phonon band is strongly dependent 
on the local temperature and can be used as a micro- 
scopic thermometer). The heat transfer coefficient be- 
tween graphene and the Si02 substrate was determined 
by modeling the heat flow using the standard heat flow 
equation with the heat transfer coefficient as the only 
unknown quantity. The authors found that using a con- 
stant (temperature independent) heat transfer coefficient 
a « 2.5 X 10^ W/m^K resulted in calculated temperature 
profiles in the graphene in good agreement with experi- 
ment. 

The heat transfer coefficient between graphene and a 
perfectly flat Si02 substrate has not been measured di- 
rectly, but measurements of the heat transfer between 
carbon nanotubes and sapphire by Maune et al[3| in- 
dicate that it may be of order a ?» 8 x 10® W/m^K. 
This value was deduced indirectly by measuring the 
breakdown voltage of carbon nanotubes, which could be 
related to the temperature increase in the nanotubes. 
Molecular dynamics calculations [4] for nanotubes on 
amorphous Si02 gives a « 3 x 10^ W/m^K, where we 
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FIG. 1: Two elastic solids with nominally flat surfaces 
squeezed together with the nominal pressure po- The heat 
current Jz(x) at the contacting interface varies strongly with 
the coordinate x = (a;, y) in the xy-plane. The average heat 
current is denoted by Jo — (Jz(x)). 



have assumed the contact width between the nanotube 
and the substrate to be 1/5 of the diameter of the nan- 
otube. A similar value can be deduced from theory [5]. 

Here we present a theoretical study of the different 
mechanisms which determine the heat transfer between 
graphene and an amorphous Si02 substrate. We study 
the contribution to the heat transfer not just from the 
area of real contact, but also the heat transfer across 
the non-contact surface area, in particular the contri- 
bution from the fluctuatingelectromagnetic field, which 
surrounds all solid objects [q], and heat transfer via the 
surrounding gas. 

Consider two elastic solids (rectangular blocks) with 
randomly rough surfaces squeezed in contact as illus- 
trated in Fig. [1] Assume that the temperature at the 
outer surfaces z = —do and z — di is kept fixed at Tq 
and Ti , respectively, with Tq > Ti. Close to the interface 
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FIG. 2: The heat transfer between two solids occur via (a) 
the electromagnetic field (e.g., photon tunneling), (b) horn 
heat diffusion or ballistic energy transfer in the surrounding 
gas (or liquid), via (c) capillary bridges and (d) via the area 
of real contact. 



the heat current will vary rapidly in space, J = J(x, z), 
where x = {x, y) denotes the lateral coordinate in the 
xy-plane. Far from the interface we will assume that 
the heat current is constant and in the z-direction, i.e., 
J = Jqz. We denote the average distance between the 
macro asperity contact regions by A (see Ref. [7[). We 
assume that X << L, where L is the linear size of the 
apparent contact between the elastic blocks. The tem- 
perature at a distance ~ A from the contacting interface 
will be approximately independent of the lateral coordi- 
nate X = (x, y) and we denote this temperature by Tq and 
r{ for z = — A and z = A, respectively. The heat current 
for \z\ >> A is independent of x and can be written as 
(to zero order in A/do and X/di): 
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where kq and ki are the heat conductivities of the two 
solid blocks. We assume that the heat transfer across the 
interface is proportional to Tq — T{ and we define the heat 
transfer coefficient a so that 
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Combining (1) and (2) gives 
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This equation is valid as long as A << L and A << 
do, di. In our application the upper solid is graphene 
which is only one atom thick. But because of the very 
high in-plane thermal conductivity of graphene compared 
to the substrate (Si02) the temperature in the graphene 
will be nearly constant, and we can define the heat trans- 
fer coefficient a using (2) in this case too. 

In Fig. [2] we show the different heat transfer pro- 
cesses we discuss below. We first study the radiative 
contribution (a) to the heat transfer [y, |8|, |9[. The heat 



flux per unit area between two black-bodies separated 
by d >> dr = ch/kBT is given by the Stefan-Boltzmann 
law 
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where To and Ti are the temperatures of solids 1 and 
2, respectively, and c the light velocity. In this limiting 
case the heat transfer between the bodies is determined 
by the propagating electromagnetic waves radiated by 
the bodies and does not depend on the separation d be- 
tween the bodies. Electromagnetic waves (or photons) 
always exist outside any body due to thermal or quan- 
tum fluctuations of the current density inside the body. 
The electromagnetic field created by the fiuctuating cur- 
rent density exists also in the form of evanescent waves, 
which are damped exponentially with the distance away 
from the surface of the body. For an isolated body, the 
evanescent waves do not give a contribution to the energy 
radiation. However, for two solids separated by d < dr, 
the heat transfer may increase by many orders of mag- 
nitude due to the evanescent electromagnetic waves-this 
is often referred to as photon tunneling. 

For short separation between two solids with flat sur- 
faces (d << dr), the heat current due to the evanescent 
electromagnetic waves is given by[6| 
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where 
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and the reflection factor 
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where e(q, w) is the dielectric function. From (4) it fol- 
lows that if i?o(q, w) and i?i(q, w) are independent of q, 
the heat current scale as 1/d^ with the separation be- 
tween the solid surfaces. 

We now apply (4) to graphene adsorbed on a nomi- 
nally flat surface of amorphous Si02. The graphene di- 
electric function was recently calculated using the mean 
fleld approximation where the electric field acting on an 
electron is the sum of the external electric field, and 
the induced field from the other electrons [lOj. In the 
calculation it was assumed that the one-particle energy 
eigenvalues are linearly related to the wavevector k via 
fsk — S7|k|, where s = ±1 indicate the conduction and 
valence bands, and where 7 ss 6.5 eVA is a band param- 
eter. In this model the Fermi wavevector kp = (ttti)^'^. 




FIG. 3: The logarithm of the heat transfer coefficient, be- 
tween graphene and an amorphous Si02 solid with a flat sur- 
face, as a function of the separation (in nanometer) d between 
the sohds. For the temperatures T = 500, 700 and 900 K. 



and the Fermi energy Ep — ^k-p, where n is the number 
of 2D carriers (electron or hole) per unit area. Note that 
in the absence of doping (and for vanishing gate voltage) 
n = 0. In the application below the excitation energy huj 
is of order the thermal energy /cbT' while the wavevector 
typically is of order 1/d, where d is the average sepa- 
ration between the graphene and the substrate surface. 
Under these conditions htu << jq which is equivalent to 
q/kp >> huj/E-p. In this limit and assuming we can ne- 
glect the influence of the temperature on the dielectric 
properties, for q < 2A;f the dielectric function takes the 
following simple form[10|: 



.1/2 



4e^ 



huj 



79 \ 11 



J_ 

2fcF 



,-, 1/2N 



In the preset case the heat transfer is associated with 
thermally excited optical (surface) phonon's in Si02 and 
electron-hole pairs in the graphene. That is, the elec- 
tric field of a thermally excited optical phonon in Si02 
excites an electron-hole pair in the graphene, leading to 
energy transfer. The excitation transfer occurs in both 
directions but if one solid is hotter than the other, there 
will be a net transfer of energy from the hotter to the 
colder solid. 

In Fig. [3]we show the logarithm of the heat transfer co- 
efhcient, between graphene and an amorphous Si02 solid 
with a flat surface, as a function of the separation (in 
nanometer) d between the solids. The results have been 
obtained from (4) for the temperatures T — 500, 700 and 
900 K. For d « 1 nm we obtain a « 3 x 10^ W/m^K. The 
heat transfer coefhcient deduced from the experiment [2] 
is aoxp ~ 2.5 X 10^ W/m^K, and we conclude that the 
field coupling gives a negligible contribution to the heat 
transfer. Note also that the photon tunneling contribu- 
tion to a depends strongly on the temperature (it in- 
creases with a factor of « 2 as the temperature increases 
from 400 K to 900 K), while the experimental data[2| 
could be fit with a temperature independent a. 

Let us now consider the contribution (b) to the heat 
transfer from the surrounding gas. Consider two solids 
with fiat surfaces separated by a distance d. Assume 
that the solids are surrounded by a gas. Let A be the 
gas mean free path. If d >> A the heat transfer between 
the solids occurs via heat diffusion in the gas. If d << A 
the heat transfer occurs by ballistic propagation of gas 
molecules from one surface to the other. In this case gas 
molecules refiected from the hotter surface will have (on 
the average) higher kinetic energy that the gas molecules 
refiected from the colder surface. This will result in heat 
transfer from the hotter to the colder surface. The heat 
transfer coefficient is approximately given bv|12| 



For q > 2/cf and Tiuj << jq the imaginary part of e vanish 
so there will be no contribution to the heat transfer from 
q > 2kp (note: the contribution from the pole of i?(q, to) 
is negligible). In the numerical calculations presented 
below we assume n = 10^'^ cm^^ giving kp — 0.056 A 
and Ep = 0.36 eV. 

The optical properties of (amorphous) silicon dioxide 
(Si02) can be described using an oscillator model|ll] 
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For air (and most other gases) at the normal atmo- 
spheric pressure and at room temperature A w 65 nm 
and Kgas ~ 0.02 W/mK. For contacting surfaces with 
surface roughness we get 
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The frequency dependent term in this expression is due 
to optical phonon's. The values for the parameters Coo, 
(a,aja,7a) and (6,0;;,, 7b) are given in Ref. [11]. 

If surface roughness occur so that the separation d 
varies with the coordinate x = {x, y) we have a sa (a(d)), 
where (..) stands for ensemble average, or average over 
the whole surface area. 



where (..) stand for ensemble average or averaging over 
the surface area. 

In the present application the surface separation is of 
order « nm so we can neglect the d-dependence in (5) and 
get a « Kgas/A « 3 X 10^ W/m^K, which is similar to 
the contribution from the electromagnetic coupling but 
much smaller than the heat transfer coefficient deduced 
from the experiment [2|- Note also that Kgas (and hence 
Qfgas) depends strongly on the temperature (it increases 
with a factor of « 2 as the temperature increases from 





FIG. 4: The contact region (black area) between two elas- 
tic solids observed at low (left) and high (right) magnifica- 
tion. The contact resistance depends mainly on the long- 
wavelength roughness, and can usually be calculated accu- 
rately from the nature of the contact observed at low magni- 
fication (left). 



300 K to 600 K), while the experimental data 2] could be 
fit with a temperature independent a. 

We now study the contribution (c) to the heat trans- 
fer from capillary bridges. If the solid walls are wet by 
water, in a humid atmosphere capillary bridges will form 
spontaneous at the interface in the vicinity of the as- 
perity contact regions. For very smooth and hydrophilic 
surfaces the fluid (in this case water) may occupy a large 
region between the surfaces and will then dominate the 
heat transfer between the solids. Similarly, contamina- 
tion layers (mainly organic molecules) which cover most 
natural surfaces may form capillary bridges between the 
contacting solids, and contribute in an important way 
to the heat transfer coefficient. The fraction of the in- 
terfacial surface area occupied by fluid bridges, and the 
separation between the solids in the fluid covered re gion , 
can be calculated using the theory developed in Ref. [13| . 
From this one can calculate the contribution to the heat 
transfer using: 

For the present system we do not expect capillary bridges 
to be important because the experiment was performed 
in dry nitrogen atmosphere. 

The study above shows that the contribution to a from 
processes (a)-(c) are much smaller than the observed heat 
transfer coefficient. Thus the heat transfer must be dom- 
inated by the only remaining process, namely heat flow 
via the area of real contact. Recent contact mechanics 
studies have shown that for elastic contact, the contact 
regions observed at atomic resolution may be just a few 
atoms wide, i.e., the diameter of the contact regions may 
be of the order of ~ 1 nm lj-[l6j . The heat transfer 



via such small junctions may be very different from the 
heat transfer through macroscopic sized contact regions, 
where the heat transfer usually is assumed to be pro- 
portional to the linear size of the contact regions (this 
is also the prediction of the macroscopic heat diffusion 
equation), rather than the contact area. In particular. 
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FIG. 5: An elastic film (a), and a semi-infinite elastic solid 
(b), in contact with a randomly rough substrate. The elastic 
modulus of the solids in (a) and (b) are such that the dis- 
placement h (in response to an external pressure po) into the 
"cavity" (with diameter A) is the same in both cases. 



if the typical phonon wavelength involved in the heat 
transfer becomes larger than the linear size of the con- 
tact regions (which will always happen at low enough 
temperature) the effective heat transfer may be strongly 
reduced. Similarly, if the phonons mean free path is 
longer than the linear size of the contact regions, bal- 
listic (phonon) energy transfer may occur which cannot 
be described by the macroscopic heat diffusion equation. 
However, as shown in Ref. [8|, for macroscopic solids the 
thermal (and electrical) contact resistance is usually very 
insensitive to the nature of the contact regions observed 
at the highest magnification, corresponding to atomistic 
(or nanoscale) length scales. In fact, the heat transfer 
is determined mainly by the nature of the contact re- 
gions observed at lower magnification where the contact 
regions appear larger (see Ref. |8| and Fig. S]). Thus, 
the thermal contact resistance of macroscopic solids usu- 
ally does not depend on whether the heat transfer occur 
by diffusive or ballistic phonon propagation, but rather 
the contact resistance is usually determined mainly by 
the nature of the contact regions observed at relative low 
magnification. 

If two semi-infinite elastic solids with randomly rough 
surfaces are squeezed in contact at the (nominal) squeez- 
ing pressure po, the heat transfer coefficient [8|] 
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The effective heat transfer coefficient k is defined by 
11 1 

where kq and ki are the heat conductivities of solid 
and solid 1, respectively. The effective elastic modulus 
E* is defined by 
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where Eq and vq are the Young's elastic modulus and 
the Poisson ratio, respectively, for solid 0, and similar 



for solid 1. The length parameter uq in (6) can be calcu- 
lated from the surface roughness power spectrum C{q) as 
described in Ref. [l7| . In general uq is of order the root- 
mean-square roughness of the (combined) surface profile. 
We have calculated uq w 2.5 nm for an amorphous Si02 
surface, with the root mean square roughness 2.5 nm 
when measured over a surface area lO/im x lO/im. 

Eq. (6) for the heat transfer via the area of real con- 
tact is valid for the case of two semi-infinite elastic solids 
in non-adhesive contact [8]. However, the present appli- 
cation involves an atomically thin elastic film in adhesive 
contact with a substrate. It may be possible to extend 
the theory presented in Ref. [8| to this case too, but here 
we will instead present an approximate treatment based 
on the theory presented above. 

In the derivation of Eq. (6) it is assumed that the 
elastic solids have a thickness much larger than the av- 
erage separation between the macroasperity contact re- 
gions (see Ref. [8]). This condition is not valid for 
graphene which is only one atom thick. Nevertheless, 
using the very simple arguments presented below, one 
can also apply the formula (6) for the heat transfer coef- 
ficient between graphene and Si02 . We first note that be- 
cause of the high (in-plane) heat conductivity of graphene 
(k w 5000 W/mK at room temperature) , the contact re- 
sistance will arise on the Si02-side of the interface, i.e., 
we can use k — ki (where ki « 1 W/mK is the heat 
conductivity of amorphous Si02) in (6). In the present 
application there is no external applied pressure but the 
graphene is bound to the Si02 substrate by adhesion. 
In a first approximation we can consider the interfacial 
interaction as the sum of the long-ranged van der Waals 
interaction and a short ranged repulsion in the contact re- 
gions, where the electron clouds of the graphene and the 
substrate overlap. We can apply (6) approximately to 
this situation if the pressure po is taken as the (average) 
force per unit aria arising from the van der Waals inter- 
action. For the separation d « 1 nm one can estimate [19| 
the van der Waals pressure to be of order po w 10^ Pa. 

Finally, we need the effective elastic modulus E* to 
be used in (6). In the present case we can neglect the 
deformations of the Si02 substrate and only include the 
deformations of the graphene. We determine E* as fol- 
lows: Assume first that the surface of a semi-infinite 
solids (with Young's modulus E*) is deformed so that 
it penetrate a distance h into a substrate cavity with di- 
ameter A, see Fig. Mjo). This cost the elastic energy 
w E*X^{h/X)'^. Here we have used that in a volume ele- 
ment ^ A'^ the typical strain is h/ X. Let us now instead 
deform the graphene so it penetrate the same distance 
h into the cavity, see Fig. [5ja). This requires stretch- 
ing the graphene 1 18| (i.e., in plane deformation) with 
i5A « X{h/Xf' over the area A^. The stored elastic energy 
is thus EXH{SX/X)^ « EX^t{h/Xf, where E sa lO^^ Pa is 
the Young's modulus of in-plane deformation of graphite 
and t is the thickness of graphene (which is equal to the 



layer spacing in graphite or i « 0.34 nm). We define E* 
so that the two elastic energies are equal which gives 
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Next note that the energy to deform the film by a dis- 
tance h can also be written as ~ pohX^, where pq is the 
applied pressure. Thus we get pohX^ « E*X^{h/X)^ or 
h/X — p/E*. Substituting this in (7) gives 
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Substituting (8) in (6) gives 
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This equation shows that a does not depend sensitively 
on Pq and A, which are not accurately known. 

Eq. (8) shows that the effective elastic modulus de- 
pends on the length scale A which is expected as the 
long-range elastic properties of the effective elastic solid 
and the graphene are different. However, we can esti- 
mate a typical A using A « (AqAi)^' ^, where Aq and Ai are 
the wavelength of the longest and shortest surface rough- 
ness components. The former is determined by the roll- 
off wavevector of the surface roughness power spectrum 
which is typically go ~ 3 x 10^ m~^ giving Aq = 27r/go ~ 
200 nm. The latter is of order a Ai « 1 nm. Thus we get 
A « 14 nm. Using (8) we obtain E* !^ 2 x 10® Pa. Using 
Po ~ 10^ Pa, E* « 10® Pa, uo ~ 3 nm and k = 1 W/mK 
in (6) gives a w 3 x 10^ W/m^K which is very close to 
the observed value. In addition, since the thermal con- 
ductivity of amorphous Si02 is only weakly dependent 
on temperature for 300 K < T < 900 K (which result 
from the short phonon mean free path in the disordered 
Si02), the heat transfer coefficient will be nearly temper- 
ature independent in the studied temperature interval, in 
agreement with the experimental results of Ref. [4I . On 
the other hand the contribution to a from the surround- 
ing gas, and from photon tunneling, depends strongly on 
the temperature. 

To summarize, we have studied theoretically the heat 
transfer between graphene and an amorphous Si02 sub- 
strate. We have found that most of the heat energy flows 
through the area of real contact, while the heat flow via 
the surrounding gas, and from photon tunneling, are both 
roughly 100 times weaker. 

Acknowledgment s 

We thank P. Avouris for drawing our attention to Ref. 
y . We thank him and M. Freitag for information related 
to the same reference. B.N. J. P. was supported by Invi- 
tation Fellowship Programs for Research in Japan from 
Japan Society of Promotion of Science (JSPS). H.U. was 
supported by the Grant-in- Aid for Scientific Research B 
(No. 21310086) from JSPS. 



[1] A.K. Geim and K.S. Novoselov, Nat. Mater. 6, 183 
(2007). 

[2] M. Freitag, M. Steiner, Y. Martin, V. Perebeinos, Z. 
Chen, J.C. Tsang and P. Avouris, NANO Letters 9, 1883 
(2009). 

[3] H. Maune, H-Y Chin and M. Bockrath, Applied Physics 
Letter 89, 013109 (2006). 

[4] Z-Y Ong and E. Pop, Phys. Rev. B81, 155408 (2010). 

[5] B.N.J. Persson and H. Ueba, in preparation. 

[6] A.L Volokitin and B.N.J. Persson, Reviews of Modern 
Physics 79, 1291 (2007). 

[7] B.N.J. Persson, Surf. Sci. Rep. 61, 201 (2006). 

[8] B.N.J. Persson, B. Lorenz and A.I. Volokitin, The Euro- 
pean Physics Journal E31, 3 (2010). 

[9] The electromagnetic field coupling contribution to heat 
transfer between a nanotube and a Si02 substrate was 
studied in: S. V. Rotkin, V. Perebeinos, A.G. Petrov 
and P. Avouris, NANO Letters 9, 1850 (2009). 



[10] E.H. Hwang and S.Das Sarma, Phys. Rev. B75, 205418 

(2007). 
[11] D.Z.A. Chen, R. Hamam, M. Soljacic, J.D. Joannopoulos 

and G. Chen, Applied Physics Letters 90, 181921 (2007). 
[12] M. Bahrami, M.M. Yanavovich and J.R. Culham, Journal 

of Thermophysics and heat transfer 18, 326 (2004). 
[13] B.N.J. Persson, J. Phys.: Condens. Matter 20, 315007 

(2008). 
[14] C. Yang, U. Tartaglino and B.N.J. Persson, Eur. Phys. J 

E19, 47 (2006). 
[15] S. Hyun, L. Pei, J.F. Mohnarie and M.O. Robbins, Phys. 

Rev. E70, 026117 (2004). 
[16] Y.F. Mo, K.T. Turner and I. Szlufarska, Nature 457, 

1116 (2009). 
[17] C. Yang and B.N.J. Persson, J. Phys.: Condens. Matter, 

20, 215214 (2008) 
[18] G. Carbone, L. Mangialardi and B.N.J. Persson, Phys. 

Rev. B70, 125407 (2004). 
[19] L.W. Bruch, Phys. Rev. B72, 033410 (2005). 



